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I.  INTRODUCTION 


The  use  of  separate  codes  for  computing  the  inviscid  flow  and  turbulent 
boundary-1 ayer  development  over  yawed  bodies  of  revolution  has  yielded  some 
very  good  solutions  for  cone  and  ogive-cylinder  shapes1.  However,  the  authors 
have  found  that  application  of  these  techniques  to  bodies  with  boattail  after¬ 
bodies  has  not  yielded  satisfactory  results  even  at  small  angle  of 
attack  (a  <  4°). 

Several  recent  publications  have  reported  supersonic  flow- field  computa¬ 
tions  using  Parabolized  Navier-Stokes  (PNS)  techniques.  These  publications 
have  reported  very  good  results  for  cone  models  for  laminar  and  turbulent 
viscous  flow2’3’4,  and  for  ogive-cylinder  bodies  and  laminar  viscous  flow4. 
The  PNS  method  appears  to  offer  an  attractive  technique  for  computing  flow 
over  bodies  with  discontinuities  in  surface  curvature  (such  as  occurs  at  the 
junction  between  the  cylinder  and  the  boattail)  since  the  inviscid  flow  and 
viscous  layer  are  computed  simultaneously.  Further,  the  PNS  method  permits 
adequate  flow-field  resolution  to  be  achieved  with  very  reasonable  computer 
costs.  This  report  describes  the  results  of  detailed  comparisons  of  PNS 

computational  results  to  experimental  measurements  for  surface  pressures  and 
turbulent  boundary- layer  profile  characteristics  of  an  ogi ve-cyl i nder-boattai 1 
body  at  Mach  =  3  and  angles  of  attack  up  to  10°.  The  PNS  code  used  is  that 
reported  by  Schiff  and  Steger5. 


II.  GOVERNING  EQUATIONS  AND  NUMERICAL  SCHEME 

A  body-conforming  £,  o,  C,  coordinate  system  (Figure  1)  is  used  which 
maps  the  body  surface  and  outer  boundary  of  the  flow  region  in  physical  space 
onto  coordinate  surfaces  of  the  computational  space.  This  transformation 


1.  Sturek,  W.  B.,  Dwyer ,  H.  A.,  Kayser,  L.  D.,  Nietubicz ,  C.  J .,  Reklis,  R. 

P.,  and  Opalka ,  K.  D. ,  "Computations  of  Magnus  Effects  for  a  Yawed, 
Spinning  Body  of  Revolution, "AIAA  Journal,  Vol  16,  No.  7,  July  1978,  pp. 
687-692.  - 

2.  Lin,  T.  C.,  and  Rubin ,  S.  G .,  "Viscous  Flow  Over  a  Cone  at  Moderate 
Incidence:  I-Hypersonic  Tip  Region, "International  Journal  of  Computers 
and  Fluids,  Vol  1,  1973,  pp.  37-57 . 

3.  Lubard,  S.  C.,  and  Helliwell,  W.  S.,  "Calculation  of  the  Flow  on  a  Cone  at 
High  Angle  of  Attack,  "AIAA  Journal,  Vol  12,  July  1974,  pp.  965-974. 

4.  Rakich,  J.  V.,  Vigneron,  Y.  C.,  and  Agarwal,  R.,  "Computation  of 
Supersonic  Viscous  Flows  Over  Ogive-Cylinders  at  Angle  of  Attack,"  AIAA 
Paper  No.  79-0131,  17th  Aerospace  Sciences  Meeting,  January  1979. 

5.  Schiff,  L.  B.,  and  Steger,  J.  L.,  "Numerical  Simulation  of  Steady 

Supersonic  Viscous  Flow,"  AIAA  Journal,  Vol.  18,  No.  12,  December  1980, 
pp.  1421-1430.  - 
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simplifies  the  application  of  surface  boundary  conditions  and  permits  the 
approximation  of  neglecting  streanwise  and  circumferential  viscous  terms  in 
hi gh-Reynol ds-number  flow  (see  Reference  5).  The  resulting  steady  thin-layer 
PNS  equations  can  be  written  in  strong  conservation-law  form  in  terms  of 
nondimensional  variables  (Reference  5)  as 


*\C  *  * 

S  IF  _3G  _1  IS 

7F  3n  3;  -  Rg  3? 


(1) 


where 


5  -  «(x) 
n  =  n(x,y,z) 
?  =  ;(x,v,z) 


is  the  axial  (marching)  coordinate 
is  the  circumferential  coordinate 
is  the  normal  coordinate 


The  inviscid  flux  vectors  in  Eq.  (1)  are 


P:J 

pV 

pW 

_1 

J  * 

pu'J^x? 

pvU 

,  F  =  J'1 

puV+nxp 

pvV+n  p 
y 

,  G  =  J"1 

puW+c  p 
PvW+CyP 

pwU 

pwV+n2P 

pwW+;zp 

(e+Ps)U 

(e+p)V 

_(e+p)W  _ 

(?) 


with  contravariant  velocity  components 


U  * 

V  =  r,  u  n  v  +  n  w 
x  y  4 

w  =  cx'J  .+  c  v  +  ;2w 


(3) 


The  internal  energy  of  the  gas  e^  is  defined  in  terms  of  the  conservative 
variables  as 


e.  =  (e/p)  -  0.5(u2  +  v2  +  w2)  (4) 

while  the  equation  of  state  for  a  perfect  gas  with  ratio  of  specific  heats  y 
i  s 


p/p  =  {y  -  l)e.  =  a2/y  (5) 

Variations  of  body  geometry  are  included  in  Eqs.  (l)-(3)  through  the  presence 

of  metric  terms  5  ,  n  ,  n  ,  etc.,  and  the  Jacobian  of  the  transformation  J, 
x  x  y 

which  appears  in  the  flux  vector. 
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The  thin-layer-model  viscous  terms,  valid  for  high-Reynolds-number  flow 
and  body-conforming  coordinates  are 


0 


u(?x+?y+?  zK  +  (^)(«A+SV?zwC)?x 
u(?x+?y+5z)v?  +  (u/3)(?xu?+?yvc+?zw?)Cy 

+  (p/3)(svsv?zw^z 

{ ( CxHyH2 ) [ (p/2 ) (u2+v2+w2) c 
+  <Pr-1 (y-1)”1 ( a2 ) ^ ]  +  (y/3)(cxu 

_  +  V  +  ?zw><?xu?  +  y5 +  ?zV} 


(6) 


For  turbulent  flow  computations  the  coefficient  of  molecular  viscosity  (p)  and 
thermal  conductivity  (k)  appearing  in  Eq.  (6)  are  computed  using  the  two-layer 
Cebeci- type  eddy  viscosity  model  reported  by  Baldwin  and  Lomax^.  The  various 
constants  within  the  model  were  set  to  the  values  suggested  in  Reference  6, 
with  the  exception  that  the  turbulent  Prandtl  number  Pr^.  was  set  to  0.8. 

Equation  (1)  is  parabolic-like  with  respect  to  £,  and  can  thus  be  marched 
downstream  in  the  5  direction  from  an  initial  data  plane  (subject  to  appro¬ 
priate  body  and  free-stream  boundary  conditions),  under  those  conditions  where 
the  local  flow  is  supersonic.  By  evaluating  the  pressure,  p$,  which  appears 

in  the  E$  flux  vector  using  the  subsonic  layer  approximation,  Eq.  (1)  can  be 

kept  stable  for  marching  for  subsonic  points  as  well.  If  p$  is  set  equal  to 

the  local  pressure  for  supersonic  points,  and  is  evaluated  from  9p  /a?  =  0 

(Figure  2)  for  points  within  the  subsonic  viscous  layer  adjacent  to  a  wall, 
Eq.  (1)  can  be  stably  marched  for  all  flows  where  U  >  0;  that  is,  for  flows 
without  streamwise  reversal  (see  Reference  5  for  associated  stability 
analysis). 

The  numerical  algorithm  used  to  advance  Eq.  (1)  downstream  in  5  is  a 
noniterative,  implicit,  approximately  factored  finite-difference  scheme. 


6.  Baldwin,  B.  S.,  and  Lomax,  H.,  "Thin  Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-25?,  16th 
Aerospace  Sciences  Meeting,  January  1978. 
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analogous  to  the  one  developed  by  Beam  and  Warming7  for  the  solution  of  tie 
unsteady  Navier-Stokes  equations.  The  algorithm  can  be  written  in  so-called 
delta  form  as 


[AsJ  +  (1  -  ct)  A<5(SnBJ)  J(AsJ)_1 
x  [A  J  +  (1  -  a)A£(6  CJ  -  Re"*  SMJ)]Aq^ 

i>  C,  Q 

=  -(A.J  -  A;j_1)qJ’  +  a(i  J  -  i^'1) 
s  s  x  s  s' 

-  (1  -  a)A5{6n[n^+1(E/J)J  +  nj+1(F/J)J  +  n^+1(G/J)J] 
+  \[^+1(E/J)j  +  ^+1(F/J)j  +  <,Jz+1(G/J)J'j 
-  rV1^}  -  [(Cx/J)J’+1Epj  -  UX/J)JV_1]  +  Pqj 

~  -1 

where  q  =  J  (p,  pu,  pv,  pw,  e) 


(?) 


The  <5's  represent  second  order  central  difference  operators  while  A 
represents  a  conventional  forward  difference.  The  Jacobian  matrices  A.  B,  and 
3  E  3  F  3G 

C  are  defined  as  and  ~  ,  respectively.  The  coefficient  matrix  M  's 

obtained  from  the  Taylor  series  linearization  of  the  viscous  vector  S.  The 
algorithm  shown  in  Eq.  (7)  is  second  order  accurate  in  Z  for  a  =  1/3,  and  is 
first  order  accurate  in  i  for  a  =  0.  The  fourth  order  dissipation  term  is 
represented  by  V  which  is  added  to  damp  high-frequency  oscillations. 

The  algorithm  is  conservative,  of  second-order  accuracy  in  the  marching 
direction,  and  can  be  either  second-  or  fourth-order  accurate  in  the  cross- 
flow  plane.  •  The  algorithm  has  been  applied  to  compute  a  variety  of  laminar 
and  turbulent  viscous  flows  and  the  results  have  been  in  excellent  agreement 
with  those  obtained  from  more  costly  time-dependent  computations.  Full 
details  of  the  notation,  of  the  PNS  assumption  and  derivation  of  the 
algorithm,  and  the  associated  stability  analysis  and  application  of  boundary 
conditions  are  found  in  Reference  5. 


7.  Beam,  R .,  and  Warming,  R.  F.,  "An  Implicit  Factored  Scheme  for  the 

Compressible  Navier-Stokes  Equations,"  AIAA  Journal,  Vol.  16,  No.  4,  April 
1978,  pp  393-402.  - 
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III.  CONICAL  INITIAL  SOLUTIONS 

In i  general,  the  initial  data  plane  for  the  marching  method  must  be  sup- 
plied  from  an  auxiliary  computation.  However,  when  treating  the  flow  over 

H°;/01nnted  bofldies>1the  aching  code  can  be  used  to  generate  its  own 
rntnal  data  plane.  As  outlined  in  Reference  5,  for  inviscid  conical  flows  a 

J!ind  1S  lecte1d  ?nd  *he  flow  variables  are  initially  set  to  free- 

anl  ^ftPr  Srh  dhe  SfuUtl0nn  1-S  marched  downstream  from  an  initial  station 

statin^  wh  h  t6Ph  the  !olu^1on  1s  scaled  to  place  it  back  at  the  original 
station.  When  no  change  in  the  flow  variables  occur  with  further  marching 

the  variables  are  constant  along  rays,  and  a  conical  solution  has  been  gene?- 

hP  w  n  the  f]ow  variables  within  the  viscous  layer  can  also  be  assumed  to 

viscous  conical  solution0,".9  rayS’  ^  ^  Pr°Cedure  «»  used  to  generate 

The  conical  viscous  approximation  is  subject  to  criticism  since  it 
assumes  the  boundary- layer  displacement  thickness  to  vary  locally  linearly 
with  distance  downstream,  instead  of  the  actual  x°*5  (laminar)  or  x0*8  * 

grownth;.  To  assess  the  accuracy  of  the  starting  procedure,  turbu- 
te  L"1' 5COiU,%so1 lut'°J1s  nwere,  generated  using  the  conical  assumption  at  two 
stations,  15.2mm  and  61.0mm  downstream  from  the  apex  of  a  10°  cone  at 

«  •  >  a  =  2.0°,  and  Re  =  2.53  x  106 */m.  The  upstream  (x  =  15.2  mm)  solu- 

52  'f  nr  a5  for  tu6  marching  code  to  obtain  a  marching  solu- 

tion  for  15.2mm  <  x  <  61.0  mm.  The  marching  solution  at  x  =  61.0  mm  free 

JnniVaT6  .errors  °f  the  conical  approximations,  are  compared  to  the  locally 
conical  viscous  flow  generated  at  x  =  61.0  mm,  in  Figures  3-4.  The  circumfer¬ 
ential  surface  pressure  distribution  is  shown  in  Figure  3,  while  boundary- 
ayer  velocity  profiles  on  the  windward  and  leeward  rays  are  shown  in  Figure 
hi  thelt  f  °W  co.ndltlons  the  conical  solution  is  virtually  indistinquish- 

thpewr0imi  the  [!iar^hln9  code  solution,  thus  demonstrating  the  suitability  of 
the  locally  conical  approximation.  y  or 


IV.  RESULTS 

A*  Model  Geometry  and  Experimental  Measurements. 

,ro  dl.menf ?ons  °l  the  ogive-cylinder-boattail  model  used  for  this  studv 

are  shown  in  Figure  5.  The  model  is  6  calibers  long  with  a  1-caliber  7° 
oattail,  and  closely  resembles  a  modern  low-drag  artillery  projectile. 

A  number  of  wind-tunnel  experiments  have  been  conducted  for  this  model 
geonetry  ,n  order  to  obtain  data  for  comparison  to  numerical  computations 
The  data  acquired  include  measurements  of  wall  static  pressure8,  turbulent 


5‘  t/"  afASt^ek>  W *  S*>  "Surface  Pressure  Measurements  on 

Slender  Bodies  at  Angle  of  Attack  in  Supersonic  Plow,"  U.S.  Army 

Ball%st%c  Research  Laboratory,  ARRADCOM  Memorandum  Report  ARBRL-MR-02876 

Aberdeen  Prov%ng  Ground,  MD,  November  1978.  AD  A064097.  * 
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boundary-layer  velocity  profiles9’10,  surface  skin  friction9,  aerodynamic 
forces,  and  flow  visualization.  The  measurements  of  primary  interest  here  are 
the  wall  static  pressure  and  boundary-layer  profile  characteristics. 

The  test  conditions  were  M  =  3  with  a  tunnel  total  pressure  of  0.298  MPa 
and  tunnel  total  temperature  of  308°K.  These  conditions  produced  a  free- 
stream  Reynolds  number  of  7.3  x  106  based  on  the  model  length.  The  boundary 
layer  was  tripped  near  the  tip  of  the  model  to  produce  a  reliable  turbulent 
flow.  All  tests  were  performed  using  SSWT  Number  One  at  the  U.S.  Army 
Ballistic  Research  Laboratory.  This  facility,  which  is  no  longer  in  opera¬ 
tion,  was  a  continuous  flow  tunnel  with  a  flexible  plate  nozzle.  The  test 
section  size  was  330  x  380mm  (13  x  15  in.). 

The  surface  pressure  measurements  were  made  using  a  model  with  a  single 
row  of  10  pressure  taps  along  the  body  length.  Data  were  obtained  as  a  func¬ 
tion  of  circumferential  position  by  rotating  the  model  in  increments  of  10°. 

Boundary-layer  velocity  profiles  were  obtained  from  measurements  of 
impact  pressure  using  a  flattened  total  head  probe.  Since  measurements  were 
not  made  of  the  total  temperature  through  the  boundary  layer,  the  temperature 
distribution  was  estimated  using  the  Crocco  linear  relationship. 

The  Preston  tube  technique  was  used  to  obtain  measured  values  for  wall 
shear  stress  on  the  cylinder  portion  of  the  model.  These  measurements  are  not 
expected  to  be  more  accurate  than  +15%  for  this  application;  however,  the  data 
are  of  qualitative  interest. 

Further  details  of  the  experiments  including  tabulations  of  the  data  are 
available  in  References  8-10. 


B.  PNS  Computations  and  Code  Performance 

Computations  were  performed  for  a  body  having  the  same  geometric  shape  as 
the  experimental  model,  and  for  flow  conditions  duplicating  that  of  the  exper¬ 
iment.  The  tip  of  the  ogive  was  replaced  with  a  cone  tangent  to  the  ogive  at 


9.  Kayser,  L.  D .,  and  Sturek ,  W.  B.,  "Experimental  Measurements  in  the 
Turbulent  Boundary  Layer  of  a  Yawed,  Spinning  Ogive-Cylinder  Body  of 
Revolution  at  Maah  3.0.  Part  II.  Data  Tabulation ",  U.S.  Army  Ballistic 
Research  Laboratory /ARRADCOM  Memorandum  Report  ARBRL-MR-02813 ,  Aberdeen 
Proving  Ground,  MD,  March  1978.  AD  A0SS707. 

10.  Kayser,  L.  D .,  and  Stv.rek,  W.  B.,  "Turbulent  Boundary  Layer  Measurements 
on  the  Boattail  Section  of  a  Yawed,  Spinning  Projectile  Shape  at  Mach 
3.0",  U.S.  Army  Ballistic  Research  Laboratory/ARRADCOM  Memorandum 
Report  ARBRL-MR-02880,  Aberdeen  Proving  Ground,  MD,  November  1978. 

AD  A065355. 


12 


x  -  15.2mm  (see  Figure  5).  Turbulent  conical  solutions  were  generated  at  that 
station  and  used  as  initial  data  for  the  PNS  marching  code. 

The  computations  presented  here  resulted  from  a  large  number  of  numerical 
runs  in  which  the  step  size,  grid  spacing,  and  smoothing  constants  were  varied 
widely.  The  step  size.  Ax,  was  varied  from  0.30mm  to  0.76mm  without  signifi¬ 
cant  changes  being  observed  in  the  solution.  The  most  critical  variable  is 
the  grid  spacing  across  the  viscous  layer.  Experience  has  shown  that  adequate 
resolution  of  the  viscous  layer  results  if  a  value  of 

■f-  _  .  +  - 

y  5  (y  =  pwuTy/V  Ut  =  /tw/pw  )  is  obtai"ed  at  the  first  point  above  the 

body  surface  at  all  axial  stations.  The  present  computations  used  a  qrid 
consisting  of  36  ci rcumferential  points  (A<f>  =  10°)  and  50  points  radially 
between  the  body  and  the  outer  boundary.  Computation  time  on  a  CDC  7600 
computer  is  2.3  sec/step  with  this  size  grid. 

C •  Comparison  Between  Computation  and  Experiment 

1#  Surface  Pressure.  The  PNS  computations  are  compared  to  experimental 
measurements,  and  to  inviscid  flow  computations  made  using  codes11’12  based  on 
MacCormack's  predictor-corrector  technique  (Figures  6-10).  Longitudinal 
surface  pressure  distributions  along  the  windward  and  leeward  rays  are  shown 

ih  FiMcreS^6"8  tor  angles  of  attack  of  4.2°,  6.3°,  and  10.4°,  respectively, 
he  PNS  and  inviscid  computations  are  both  in  excellent  agreement  with  the 
experimental  measurements  for  angles  of  attack  up  to  4.2°.  However  it  is  of 
interest  to  note  that  the  PNS  computations  exhibit  consistently  better  agree¬ 
ment  with  experiment  in  the  vicinity  of  the  discontinuities  in  streamwise 
surface  curvature  at  the  ogive-cylinder  and  cyl inder-boattai 1  junctions. 

.  “  -  b*i  1  .gure  7)  tbe  pNS  computation  continues  to  exhibit  good  agreement 
with  the  experimental  measurements,  while  the  inviscid  computation  indicates  a 
discrepancy  on  the  leeward  side  of  the  boattail.  The  discrepancy  is  attri¬ 
buted  to  the  strong  viscous  interaction  caused  by  the  appearance  of  crossflow 
separation  and  the  roll  up  of  leeward-side  vortices  in  the  vicinity  of  the 
boattail.  As  the  angle  of  attack  is  increased  further,  the  crossflow  separa- 
I°"  p®910"  rextends  Progressively  farther  forward  along  the  body.  At  a  = 
1U.4  the  PNS  computation  exhibits  good  agreement  with  the  experimental  mea¬ 
surements,  while  the  inviscid  results  deviate  from  the  experiment  on  both  the 
leeward  side  of  the  cylinder  and  the  boattail. 

The  appearance  of  a  crossflow  separation  region  is  further  indicated  by 
the  circumferentia1  surface  pressure  distributions.  Examples  of  comparisons 
of  these  distributions  are  shown  for  a  =  6.3°  in  Figure  9  and  for  a  =  10.4°  in 
Figure  10  at  two  longitudinal  stations;  one  on  the  cylinder  portion  of  the 


11.  Schiff,  L.  B.,  "Nonlinear  Aerodynamics  of  Bodies  in  Coning  Motion ,"  AIM 
J.oum.al,  Vol .  10,  No.  11,  November  1972,  pp.  1517-1522 . 

12.  Sanders,  B.  R.,  and  Dwyer,  H.  A.,  "Magnus  Forces  on  Spinning  Supersonic 
Cones,  Part  II:  The  Inviscid  Flow, »  A  UA  IquxmlL  Vol.  14  No.  5 

May  1976,  pp.  576-582.  * 
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model  near  the  boattail,  the  second,  midway  on  the  boattail.  At  a  =  4.2°  such 
comparisons  (not  shown)  indicated  uniformly  good  agreement  between  both  compu¬ 
tations  and  the  experiment  for  all  body  stations.  At  a  *  6.3°  (Figure  9a)  the 
comparison  on  the  cylinder  indicates  excellent  agreement  between  the  PNS 
computation  and  experiment  and  the  appearance  of  a  systematic  discrepancy 
between  the  inviscid  computation  and  experiment  for  100°  <  <p  <  150°.  This 
trend  is  accentuated  for  flow  on  the  boattail  (Figure  9b).  The  comparison 
shown  in  Figure  in  for  a  =  10.4°  indicates  further  development  of  a  crossflow 
separation  for  flow  over  the  boattail  and  cylinder.  The  abrupt  rise  in  exper¬ 
imental  surface  pressure  at  <j>  =90°  indicates  the  location  of  the  crossflow 
0n  P°"|n^*  The  inviscid  computation  predicts  a  crossflow  shock  at 
L“  140  ^hlch  ls.  not  Present  in  the  experimental  data.  At  this  incidence  the 
KN3  computation  is  in  only  fair  agreement  with  the  experiment  and  suggests  an 
upper  limit  of  applicability  of  the  present  computational  technique  of  a  =  6° 
for  this  class  of  body  shapes. 

.  Skir  Friction.  An  example  of  results  obtained  for  the  skin-friction 
coefficients  are  shown  in  Figure  11  for  a  fixed  longitudinal  station  on  the 
cylinder  near  the  boattail,  at  circumferential  stations  around  the  body, 
omputed  results  for  a  =  4.2°  are  shown  and  are  compared  with  Preston  tube 
measurements  made  at  that  incidence.  The  computed  skin-friction  values  for 
a  =  4.2  are  less  than  the  experimental  ones  on  the  windward  side  and  greater 
than  experiment  on  the  leeward  side  of  the  body;  however,  the  comparison  is 
within  the  expected  experimental  uncertainty  of  +15%  (the  error  bars  in  Figure 
11  indicate  +0.001  cf  which  is  approximately  +10%).  In  addition,  computed 

results  for  a l-  6 .3°  are  shown  in  Figure  11  which  illustrate  a  shift  in  trend 
c  ,oth5.  com.Pu.ted  skin-friction  coefficient  with  increasing  incidence.  At  a  = 
6.3  the  minimum  value  of  cf  is  reached  at  4,  =  150°  as  opposed  to  <p  =  180°  as 

in  the  0-4.2°  case.  This  is  a  further  indication  of  the  developing  leeward- 

press  u  re" Vis  tr i but^o n sT  WUh  increasing  incidence  previously  shown  in  the 

3’  Jtreamw1se  Velocity  Profiles.  A  more  sensitive  test  of  the  accuracy 
of  the  computational  technique  applied  to  this  flow  is  the  comparison  of 
measured  and  computed  boundary- layer  velocity  profiles.  Such  comparisons  are 
shown  in  Figures  12-16  for  two  longitudinal  stations;  station  A  on  the 
cylinder  near  the  boattail,  and  station  B  on  the  boattail  (see  Figure  5). 
Each  figure  shows  the  velocity  profiles  for  a  particular  longitudinal  station, 
at  circumferential  stations  ranging  from  the  windward  to  leeward  ray  in  30° 
increments.  The  nondimensi onal  streamwise  velocity  components,  u,  are  plotted 
versus  physical  distance  y  measured  radially  from  the  body  surface  in  milli¬ 
meters,  rather  than  against  normalized  y/6.  This  method  of  plotting  prevents 

scaling  differences  between  the  computation  and  experiment  from  givinq  a  false 
comparison.  a  y  c 


Comparisons  for  o  -  4.2°  are  shown  in  Figure  12  for  the  cylinder  (Station 
and  in  Figure  13  for  the  boattail  (Station  B).  The  agreement 
general,  excellent.  This  is  particularly  true  on  the  windward  side  . 
body.  However,  a  slight  discrepancy  is  visible  for  the  profiles  at  d>  =  150° 
is  more  apparent  in  Figure  13. 
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Similar  comparisons  for  a  =  6.3°  are  shown  in  Figures  14-15.  At  this 
angle  of  attack  the  windward  side  measured  and  computed  profiles  are  still  in 
excellent  agreement.  However,  the  discrepancy  between  the  profiles  at  <|>  = 
150°  is  substantial,  particularly  at  the  boattail  station  (Figure  15). 

Note  that  this  discrepancy  is  less  strongly  reflected  in  the  surface  pressure 
distribution  at  the  corresponding  station  (Figure  9b). 

To  assess  the  relative  roles  of  ci rcumferential  pressure  gradients  versus 
the  expansion  over  the  boattail  as  the  source  of  the  discrepancy  in  the  veloc¬ 
ity  profiles,  a  comparison  was  also  made  at  a  =  6.3°  for  a  body  with  a 
straight  cylindrical  afterbody  replacing  the  boattail.  This  comparison  is 
shown  in  Figure  16  at  a  longitudinal  station  comparable  to  that  of  the  boat- 
tail.  The  same  computation  overprediction  of  streamwise  velocity  at  <{>  =150° 
is  present  as  well  observed  for  the  boattailed  body.  This  result  suggests 
that  the  source  of  the  discrepancy  is  not  flow  expansion  over  the  boattail, 
but  rather  is  linked  primarily  to  the  development  of  the  leeward-side  vortex 
structure. 

4.  Ci rcumferential  Velocity.  A  descriptive  indication  of  the  develop¬ 
ment  of  the  leeward-side  vortex  structure  with  increasing  incidence  can  be 
seen  in  Figures  17-19.  These  figures  show  the  projections  of  the  computed 
flow  velocity  vectors  onto  the  crossflow  plane  (i.e.,  a  plane  normal  to  the 
body  axis  of  symmetry)  located  at  x  =  324mm,  on  the  boattail,  for  angles  of 
attack  of  4.2°,  6.3°,  and  10.4°,  respectively.  At  a  =  4.2°  no  crossflow 

separation  (characterized  by  ci rcumferential  flow  from  the  leeward  toward  the 

windward  side  of  the  the  body)  is  observed  at  this  longitudinal  station. 
At  a  -  6.3°  crossflow  separation  is  observed  at  a  crossflow  separation  angle, 
<j>s,  of  132°.  This  angle,  obtained  from  interpolation  of  the  computed  circum¬ 
ferential  velocity  just  above  the  body  surface,  is  shown  in  Figure  18. 
Examination  of  additional  longitudinal  stations  indicates  that,  at  this  angle 
of  attack,  the  crossflow  separation  region  starts  just  downstream  of  the 
cyl  i nder-boattai  1  junction.  At  a  =  10.4°  the  crossflow  separation  angle  <j>  = 

117°  (Figure  19),  and  the  region  of  crossflow  separation  starts  farther  up  the 
body,  in  the  vicinity  of  the  ogi ve-cyl i nder  junction. 

Experimental  vapor-screen  flow  visualization  photos13  of  the  leeward-side 

'L°9oiC?%on  the  boat£ail  are  shown  in  Figures  20-22,  at  angles  of  attack  of 
4.2  ,  6.3  ,  and  10.4°.  In  the  vapor-screen  technique  water  vapor  is  intro¬ 
duced  into  the  wind  tunnel  to  produce  a  fog  in  the  test  section.  A  thin, 

intense  plane  of  light,  oriented  normal  to  the  body  axis  of  symmetry  at  a 
longitudinal  station  midway  along  the  boattail,  illuminates  the  fog.  The 

leeward-side  vortices  are  visible  as  dark  regions  within  the  light  plane.  The 
growth  of  the  vortices  with  increasing  angle  of  attack  is  clearly  visible. 

Comparison  of  the  vapor-screen  photos  and  the  computed  ci rcumferential 
velocity  plots  suggests  that  the  computed  regions  of  crossflow  separation  are 


13,  Uietubiez,  C .  Unpublished  BEL  Wind  Tunnel  Data  (pvivate 
aormunioation) .  1979. 
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Dlanp1imuhno\nh,-ecXtent  ^an  the  re?ions  °f  vorticlty  observed  in  the  light 


plane.  While  this  comparison  is  only  qualitative, 
previous  comparisons  of  surface  pressure  and  of 
velocity  profiles.  Although  the  present  numerical 
diction  of  crossflow  separation  it  is  obvious  that 
between  the  computation  and  experiment  develop  at 
than  6°.  The  reason  or  reasons  for  this  limitation 


it  is  consistent  with  the 
streamwise  boundary-layer 
technique  allows  the  pre- 
substantial  discrepancies 
angles  of  attack  greater 
are  not  known  for  certain; 


potential  sources,  include  the  neglect  of  circumferential  viscous  terms  within 
the  thin-layer  viscous  model,  and  uncertainties  in  three-dimensional  turbu- 
lence  modeling.  However,  it  is  of  interest  to  note  that  McRae  and  Hussaini14 
who  used  a  numerical  technique  which  retains  the  circumferential  viscous  terms 
and  employs  a  similar  eddy-viscosity  model  in  a  study  of  turbulent  cone  flows 
at  large  incidence,  also  found  signifcant  discrepancies  within  the  crossflow 

?nSathtl0wwre91-0n*  -In  their  study’  the  discrepancies  were  reduced  by  modify- 
lng^the  eddy- viscosity  turbulence  model  in  the  separated  region.  Rakich,  n 
al.  reported  similar  behavior  for  flow  over  cones. 


V.  CONCLUDING  REMARKS 

Jhius  [®P°rt  has  described  the  results  of  a  numerical  computational  study 
Parabolized  Navier-Stokes  marching  code  recently  developed  by 
Schiff  and  Steger  has  been  exercised  for  a  6-caliber,  ogive-cylinder-boatta-'l 
shape  at  incidence.  Extensive,  detailed  comparisons  to  experimental  data  at  M 

-  j,  a  <  10.4  have  been  performed  to  evaluate  the  accuracy  and  stability  of 
the  numerical  technique. 

Comparisons  have  been  made  showing  excellent  agreement  between  computa¬ 
tion  and  experiment  for  surface  pressures  at  angles  of  attack  up  to  6.3°  fit 
an  angle  of  attack  of  10.4°,  the  comparison  is  degraded  due  to’the  lability 
to  accurately  model  the  signifcant  streamwise  vortical  flow.  Excellent  agree- 

/  ?0,oachluVed  f0r  velocity  profiles  at  a  =  4.2°.  At  an  angle  of 
attack  of  6.3  the  agreement  is  good  with  the  exception  of  the  d>  =  15C° 

Ih  fat!°n*.  A?ai’n>  this  is  attributed  to  the  inability  to  accur¬ 
ately  model  the  developing  leeward-side  vortex  flow.  The  previse  cause  of  the 
inaccuracy  has  not  been  determined;  however,  three-dimensional  turbulence 

fnd  finc.1.usl0.n  of  circumferential  viscous  terms  have  been  identified 
as  areas  for  further  investigation. 

tn  Altbough  this  computation  technique  appears,  as  a  result  of  this  study, 

J  cinnVfTl6*  t-°  “  <  6  f°r. thls  class  of  bodies,  the  computed  results  do  show 
nJ  J  l 'ii  imPr0vement  ™  capability  compared  to  techniques  which  compute 

inviscid  flow  and  viscous  effects  separately  (see  Reference  1).  This  increase 
in  capability  is  achieved  without  any  significant  increase  in  computation  cost 

:uhx,cDmr7t6ofoorc„c:Prn:r  resources  9reater  than  those  current*  aM,'ab'2 
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INITIAL  DATA  PLANE 


Figure  1.  Coordinates  and  Notation 


19 


Figure  2.  Subsonic  Layer  Approximation 
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Figure  3.  Circumferential  Surface  Pressure  Distribution  on  10°  Cone 
M  =  3.0,  ct  =  2.0  ,  ReY  =  1.54  x  106 
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Figure  4.  Windward  and  Leeward  Boundary-Layer  Velocity  Profiles  on  10° 
Cone;  M  =  3.0,  a  =  2.0  ,  Rex  =  1.54  x  106 
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NOTE  .’DIMENSIONS  ARE  IN  MILLIMETRES 

A,B  BOUNDARY  LAYER  SURVEY  STATIONS 
C  CONICAL  STARTING  SOLUTION  STATION 


Figure  5.  Ogive-Cyl inder-Boattail  Model  Dimensions 
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Figure  6. 


Axial  Surface  Pressure  Distribute 
Boattail  Body;  M  =  3.0,  a  =  4.2°, 


on  on  Ogive-Cyl inder- 
Re,  =  2.13  x  107/m 
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Figure  7.  Axial  Surface  Pressure  Distribution  on  Ogive-Cyl inder- 
Boattail  Body;  M  =  3.0,  a  =  6.3°,  Re<jo  =  2.13  x  l07/m 
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Figure  8.  Axial  Surface  Pressure  Distribution  on  Ogive-Cyl inder- 
Boattail  Body;  M  =  3.0,  a  =  10.4°,  Re  =  2.13  x  107/m 

oo  ' 
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Figure  9.  Circumferential  Surface  Pressure  Distributions; 
M  =  3.0,  a  =  6.3°,  Re^  =  2.13  x  106/m 

a.  x  =  279mm  (0.915  ft),  Cylinder 
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Figure  9.  Circumferential  Surface  Pressure  Distributions; 

M  =  3.0,  a  =  6.3°,  Re  =  2.13  x  106/m 

00 

b.  x  =  330mm  (1.082  ft),  Boattail 


28 


Figure  10.  Circumferential  Surface  Pressure  Distributions; 
M  =  3.0,  a  =  10.4°,  Re^  =  2.13  x  io6/m 

a*  x  =  279mm  (0.915),  Cylinder 
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Figure  10.  Circumferential  Surface  Pressure  Distributions; 
M  =  3.0,  a  =  10.4°,  Re  =  2.13  *  106/m 

oo 

b.  x  =  330mm  (1.082  ft),  Boattail 
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Figure  11.  Circumferential  Skin  Friction  Coefficient  Distribution; 
M  =  3.0,  x  =  254mm,  Rex  =  5.41  x  106 
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0  EXPERIMENT,  REF  9 
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Figure  12. 


Boundary-Layer  Velocity  Profiles 
Boattail  Body;  M  =  3.0,  a  =  4.2°, 
x  =  254mm  (0.833  ft).  Cylinder 


on  Ogive-Cyl inder- 
Re^  =  2.13  x  l07/m 
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Figure  13.  Boundary-Layer  Velocity  Profiles  on  Ogive-Cylinder- 
Boattail  Body;  M  =  3.0,  a  =  4.2°,  Re  =  2.13  *  107/m, 
x  =  324mm  (1.063  ft),  Boattail 
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Figure  14.  Boundary-Layer  Velocity  Profiles  on  Ogive-Cyl inder- 
Boattail  Body;  M  =  3.0,  a  =  6.3°  Re  =  2.13  x  l07/m, 
x  =  254mm  (0.833  ft).  Cylinder 
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Figure  15.  Boundary- Layer  Velocity  Profiles  on  Ogive-Cyl inder- 
Boattail  Body;  M  =  3.0,  a  =  6.3°,  Re  =  2.13  x  107/m, 
x  =  324mm  (1.063  ft),  Boattail 
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Figure  16.  Boundary-Layer  Velocity  Profiles  on  Ogive-Cylinder  Body; 

M  =  3.0,  a  =  6.3°,  Re^  =  2.13  x  107/m,  x  =  318mm  (1.043  ft) 
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Figure  17.  Crossflow  Plane  Velocity  Vectors  on  Boattail;  M  =  3  0, 
“  =  4.2°,  Re^  =  2.13  x  l07/m,  x  =  324mm  (1.063  ft) 
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Figure  18.  Crossflow  Plane  Velocity  Vectors  on  Boattail;  M  =  3.0, 
a  =  6.3°,  Rem  =  2.13  x  107/m,  x  *  324mm  (1.063  ft) 
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Figure  19.  Crossflow  Plane  Velocity  Vectors  on  Boattail;  M  =  3.0, 
a  =  10.4°,  Re  =  2.13  x  107/m,  x  =  324mm  (1.063  ft) 

oo 
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Figure  20.  Vapor-Screen  Flow  Visualization  on  Boattail;  M  =  3.0, 
a  =  4.2°,  Re  =  2.13  x  107/m 

oo 
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Figure  21.  Vapor-Screen  Flow  Visualization  on  Boattail;  M  =  3.0, 
a  =  6.3°,  Re  =  2.13  *  107/m 

oo 
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Figure  22.  Vapor-Screen  Flow  Visualization  on  Boattail;  M  =  3.0, 
a  =  10.4°,  Re  =  2.13  x  107/m 

oo 
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list  of  symbols 


a 

speed  of  sound 

cf 

skin  friction  coefficient 

CP 

specific  heat  at  constant  pressure 

e 

■4 

t0by'p®"e29y  Per  unft  vo,ume  of  fluid;  normalized 

ei 

^  A  A  A  A 

Es»P»G,q 

internal  energy,  normalized  by  aw2 

flux  vectors  of  transformed  gasdynamic  equation  (Eq.  2) 

J 

Jacobian  of  transformation  between  physical  and 
computational  coordinates  P  y  d 

L 

reference  length 

M 

Mach  number 

P 

pressure,  normalized  by  p  a  2 

OO  00 

Pr 

Prandtl  number,  \i  c  /< 

“  p  oo 

Re 

Reynolds  number, 

A. 

Re 

Reynolds  number  (Eq.  1),  p  a  L/p 

Ak, 

S 

viscous  flux  vector  (Eq.  6) 

u.v.w 

Ca:Lesia?.ve!ocit^  components  along  the  x  v  z  axec 
respectively,  normalized  by  a  ’  y’  1  axes» 

J  00 

u,v,w 

Contravariant  velocity  components  (Eq.  3) 

x.y.z 

physical  Cartesian  coordinate  axes  (Figure  1) 

a 

angle  of  attack 

Y 

ratio  of  specific  heats 

IC 

COs"eaieJalSe  conductfvit>'>  normalized  by  free- 

OO 

U 

^vlfueT1  °f  ViSC0Sity’  normalized  by  free-stream 

oo 

5,n,c 

43 


LIST  OF  SYMBOLS  (Continued) 

p  density,  normalized  by  free- stream  density 

<i>  circumferential  angle  (Figure  1) 

Subscri pts 

00  free- stream  conditions 

w  body  surface  values 

x  based  on  axial  distance  from  nose 
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